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ABSTRACT 

A new summability method was tested to calculate the critical exponent v of the local- 
ization length for the symplectic case derived from the non-linear cr-model. Although we 
used the same series as Hikami and others, unlike them we were able to resum the series 
in two-dimensions (2D) and obtain the result v ~ 1. Values of v in 2 + e dimensions seem 
to saturate the Harris inequality up to e = 0.2. 
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1. Introduction - Anderson localization is known as a problem where the wave function 
localizes due to a random potential scattering [|J . The critical behaviour of the localization 
transition is described by the non-linear (NL) cr-model [||. The NLa-model explains 
succesfully various aspects of the Anderson localization, including the nonsingular density 
of states and three different universality classes - the orthogonal, unitary, and symplectic, 
depending on whether time-reversal symmetry is preserved or spin-flip scattering occurs. 

The critical exponent v describes the behaviour of the localization length A near the 
mobility edge E c , 

\~{E c -E)-». (1) 

To obtain critical exponents at the transition point in the framework of the NLu-model, 
the scheme of the minimal subtraction by the dimensional regularization is usually used. 
The critical exponent v of the correlation length or the localization length is then given 
by the relation 

where t c is the zero of the f3 function and j3' is t he derivative of (3 with respect to t [[|. 
In general, the f3 function of the NLcr-model in 2 + e dimensions is given as 

oo 

P{t)=et- y £a n t n 1 (3) 

n=2 

where the coefficients a n 's do not depend on e and a 2 ^ Q. In the symplectic case, 
which is the universality class corresponding to time-reversal symmetry and strong spin- 
orbit coupling, the f3 function to five-loop order is given by 

(3 so (t) =et + t 2 - |C(3) t 5 - ^C(4) t 6 + 0(f), (4) 

where ( is the Riemann ^-function 

2. Methods - To extract the critical exponent u, only the Borel-Pade method has been 
used so far The Borel summability method, when applied to the series of the form 
(H), consists in the following transformation, 

£a n r^~ / e-^^^du. (5) 

n=l t J ° n=l n - 

For finite series, relation (|5|) is useless since it is an identity. One can integrate term by 
term, by using 

1 roo 

- / e~ u/t u n du = n\ t n . 
t Jo 



However, for infinite series the second form in (|5]) may have a much wider range of ap- 
plicability. For example, in the case of the series zU the second form converges in 
the whole complex halfplane Re z < 1 despite that the original series being divergent for 
\z\ > 1, outside its radius of convergence. Having only a finite number of a series up to 
order N at one's disposal, such as in the present case, the integrand in (^j) can be approx- 
imated by the Pade approximation which generates an infinite power series coinciding up 
to order TV with the original series. The resulting method is called the Borel-Pade method. 

Several years ago we have developed a method || which consists in a similar transfor- 
mation as @ but with e _U// * replaced by e~ eU/t and n! replaced by 

POO 

fiOn) = / e- eU u n du. (6) 
Jo 

Moments fi(n) are increasing more slowly then n! (they behave roughly as ln n (n) when 
n — > oo), but as far as analytic properties are concerned this results in a wider region 
of convergence. If f(z) denotes an analytic continuation of a power series with a non- 
zero radius of convergence, then our method gives a finite result f(z) in the so-called 
Mittag-Leffler star of f(z) ||. In the case of the series Y^=o z n it means that our method 
converges in the whole complex plane except for the interval [1, oo). For a general power 
series, the Mittag-Leffler star is obtained in the following way (see Fig. [I]). First, one 
draws rays from the origin and passing through singularities of f(z). The Mittage-Lefffer 
star is the region which remains in the complex plane after the part of the ray beyond 
each singularity is removed. We recall that the actual region of convergence for the Borel 
method is a polygon which is obtained by removing half-planes from the complex plane 
which lie behind a perpendicular to the ray from the origin passing through the singularity 
(see Fig. [I]). If one has only a finite number of terms of a series at disposal, one can use 
the Pade approximation of the integrand in (J^) in the same way as when the Borel method 
is used. We shall call the resulting method the /z-Pade method. Both, the Borel method 
and our method are so-called analytic moment- const ant summability methods [||, §]. They 
can be applied to both convergent and divergent series. 

3. Results - One of the motivations in using the /z-Pade method for the symplectic 
case was the fact the series (f|) is not Borel summable and the Borel-Pade method does 
not work for this case M. In the 2D symplectic case, to leading order, so-called weak anti- 
localization occurs M and until recently it appeared that this result remain unchanged 
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a) 




b) 

Figure 1: a) The actual region of convergence for the Borel method in the case of the 
Taylor series of an analytic function f(z) is a polygon which is obtained by removing half- 
planes from the complex plane which lie behind a perpendicular to the ray from the origin 
passing through the singularity (a cross), b) In the case of our method, the actual region 
of convergence is obtained by only removing the part of the ray behind each singularity. 
Such region is called the Mittage-Lemer star. 
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Figure 2: Typical behaviour of the (3 function for the symplectic case after it has been 
resummed by the diagonal /z-Pade method. 

by higher order terms. This included the strange result that there is no fixed point (and 
hence no transition and no localization) for the 2D symplectic case ||. We have applied 
the diagonal [3/3] /x-Pade method directly to the f3 so {t) and looked for its zero t c (e) as a 
function of e (see Fig. g). By using relation (g) we have found that as e [ 0, the critical 
exponents 

v -> 0.98 (7) 

in the 2D symplectic case. The j3 function for the orthogonal universality class is related 
to (3 function for the symplectic universality class by substituting t in Eq. (|j) by — 2t [|J. 
Therefore, because the diagonal [3/3] /i-Pade method works for the symplectic case, it 
cannot be applied to the orthogonal case, since the Pade approximant, having a polynomial 
of order 3 in the denominator, develops a pole in the integration interval. A similar 
statement applies to the diagonal [3/3] Borel-Pade method, which, in contrast, works for 
the orthogonal case and therefore does not for the symplectic case [Q. 

In the following, we have scanned the e-expansion by varying e = d — 2, where d is the 
space dimension, within the interval e G (0, 1]. Results are presented in Table I. 
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TABLE I. (3'(t c ) and the critical coefficient v as a function of e 



e 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 




-1.05 


-1.14 


-1.27 


-1.47 


-1.74 


-2.1 


-2.58 


-3.2 




-0.50 


V 


0.95 


0.88 


0.79 


0.68 


0.57 


0.48 


0.39 


0.31 
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4. Discussion - We have checked whether our results satisfy the Harris inequality |J, 

, > §, (8) 

derived under the assumption of the validity of one-parameter scaling. Up to an error due 
to the finite number of terms for (3 so {t) [see (§)], our result seems to saturate the inequality 
AD up to e = 0.2. As £ increases further, v ceases to satisfy (§) and eventually around 
£ = 0.9 the diagonal /x-Pade method collapses. The reason is that the Pade approximant 
develops a pole on the integration interval. For £ = 1, the /i-Pade method starts to work 
again and our result 

v = 2 (9) 

satisfies inequality (H). However, because of the collapse of the diagonal /i-Pade method 
at £ = 0.9, this result must be taken with some reservations. Since the e expansion (D) 
is an asymptotic expansion, it is difficult to make an extrapolation too far away from the 
limiting point (D = 2 in our case), given that the extrapolation is based on incomplete 
knowledge of the (3 function. Nevertheless, it provides a substational improvement over 
previous analytical results, and there is always a chance that further terms of the (3 func- 
tion will make our results better. It is worthwhile to mention that the actual convergence 
or divergence of an asymptotic series is not as important for applicability of an analytic 
summability method as whether asymptotic series obey the strong asymptotic conditions 
(SAC) 0, [IIJ. The latter ensure that there exists only one function with the required 
analytic properties and a given asymptotic expansion. For example, given a convergent 
asymptotic series X)o° z ™ i n ^ ne right complex half plane, without the validity of the SAC 
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the sum of this series can be any function of the form 1/(1 — z) + As B e~ c / z with A, B, 
and C > arbitrary constants. 

In the 2D symplectic case, our result (0) for the critical exponent v is smaller than 
v ~ 2.5 Jll or z/ ~ 2.74 [12| obtained by numerical scaling analysis using, respectively, 



the Evangelou-Ziman model |13| or Ando's model ||14|| . A similar disagreement between 
field theory predictions and tight-binding scaling methods is also known to exist for the 
3D orthogonal case, where the former yields v ~ 1 and the latter v ~ 1.4. It is interesting 
to note that a disagreement also exists for the results for the critical exponent v obtained 
by the numerical scaling analysis and that obtained by the best fit to the critical level 
distribution P(s) in the 2D symplectic case at the mobility edge \ TE \, 



P{s) = J Bs 4 exp(-As 7 ). 



Here parameter 7 is given by 



7=1-7-7 



1 



(10) 



11) 



d is the dimensionality of the system, A is a numerical factor which depends on the 
dimensionality, and B is to be found from the normalization conditions j0|. The critical 
exponent u* in (|ll|) should be identical to v. However recent numerical analysis implies 
that v* = 0.7 ±0.08 pf or v* ~ 0.77 for Ando's model [0, and v* = 0.83 ±0.7 |l| 



for the Evangelou-Ziman model |13|]. A comparison with (j8|) shows that v* even violates 
the Harris inequality [however v* satisfies a weaker relation, v* > 1/d, which suffices 
to derive flTOD I. Surprisingly enough, our result (|^) for the critical exponent v in the 2D 
symplectic case for the localization length is actually very close to v* and it seems to 
be tempting to say that the critical exponent obtained from the NLcx is just v* . That, 
however, would be premature, because in the 3D orthogonal case the overall behaviour 
of P(s) has been claimed to be well fitted with v obtained from the numerical scaling 
analysis fli~8 |. 



To conclude, we have demonstrated that the /z-Pade method can be useful in deter- 
mining critical exponents. In principle, our method can be used as an alternative to 
the Borel-Pade summability method, whenever a result is obtained in a form of a power 
(asymptotic) series, such as in the case of high temperature expansion, weak coupling 
expansion, etc. In connection with disordered systems, it would be interesting to use 
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our method for a calculation of the density of states and the diffusion constant for the 
problem of an electron moving in two dimensions in the lowest Landau level and in a ran- 
dom potential, which has been analyzed by the Borel-Pade method [|19| . It was suggested 



that our summability method be used for the location of critical points ||10|| , since it finds 
singularities of an analytic function much more precisely than the Borel method (cf. Fig. 
[1]). An application of the /z-Pade method to other problems will be given elsewhere. 
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